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Abstract 

The dependence of DNA assembhes conductance on relative humidity is investigated theoreti- 
cahy. Fohowing earlier suggestions we consider the ionic conductivity through the layers of water 
adsorbed by DNA molecules. The increase of humidity results in a growing water layer. The bind- 
ing energy of ions depends on the thickness of the water layer due to change in water polarization. 
This dependence is very strong at smaller thicknesses of water layers due to the low-dimensional 
confinement of an electric field in water. We show that the associated change in ion concentration 
can explain the 6 orders of magnitude increase of conductivity, with relative humidity growing from 
0.05 to 0.95. 

PACS numbers: 
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I. INTRODUCTION 



Since the first experimental demonstration of cliarge transfer [l| 
sidered promising candidates for the reahzation of nanodevices [2 



DNA molecules are con- 
Extensive 



3l,Q,B,y,3. 



investigation of molecular junctions made of long DNA molecules has resulted in many 
interesting, but sometimes controversial observations. One common feature discovered in 

different experimental setups with long DNA molecules , Q, Q, Q, Q, Q is the dra- 
matic sensitivity of dry DNA conductivity to the humidity of air. Indeed, the increase of 
humidity from 5% to 95% raises the conductivity of DNA molecular assemblies by six orders 
of magnitude. 

The nature of this strong conductance dependence of relative humidity has been discussed 
by various authors [s], o], [lO, 11, H, Q- It was suggested sl to interpret observations 
in terms of ionic conduction through the water layers around DNA molecules. According 
to Refs. j^, Q, Q], ions HsO"^ and 0H~, formed in water self- ionization processes, are 
responsible for observed conductivity. The binding energy of water depends on the relative 
humidity because the adsorption of water increases water solubility. Alternatively, there are 
theoretical models suggesting major electronic mechanism of conductivity (see e. g. Refs. 



15 



16j). In our opinion, ions are the instruments of conduction rather than electrons. 



According to Ref. [17| the activation energy of dry DNA for the carrier hopping through 
DNA duplex is smaller than that of k- DNA in buffer, which contradicts with the increase of 



conductance with increasing relative humidity. Also according to experiments 12|, |l3| single 



and double stranded DNA do not show a dramatic difference in conductance. In case of 
electronic mechanism of conductivity this result conflicts with the absence of charge transfer 
in a single stranded DNA 1181. It is important to note that DNA conductance behaves 

n 

differently in short DNA junctions [19|, where the conductivity is most likely associated 
with the electron transport. 

In this manuscript we develop an electrostatic model to investigate the effect of humidity 
on ionic conductivity of DNA assemblies. We treat each DNA molecule independently and 
assume each is surrounded by a layer of water Fig. [T] In Sec. [TTl the humidity dependent 



thickness of the water layer is estimated using experimental data [20|, |21|, 23| and each 
dielectric layer in our model is characterized by its bulk dielectric constant. Since dielectric 
constant of water {ew = 80) is much larger than that of the air (eo = 1) the electrostatic 
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FIG. 1: Field confinement in the water tube formed around DNA in a humid environment (see 
text for details). 



interaction in the thin water layer is different from that in the bulk water. The boundary 
conditions for the electric field between two dielectrics (water and air) results in the field 
confinement inside the water layer strongly affecting ion binding energies. Field confinement 
means that the electric field of interacting ions in water is substantially trapped in the water 



layer, and this effect causes DNA conductivity to be dramatically sensitive to 



lumi dity 



changes. The field confinement has been studied earlier for ionic channels [2j, |25|, |26|, |27| , 
where we borrowed this idea. The field confinement results in stronger attraction between 
positive and negative ions. Accordingly, fewer water molecules between ions results in an 
increase in binding energy and a decrease in conductance. In Section IIIII we calculate the 
binding energy associated with electric field confinement using the exact solution of Poisson 
equations in the cylindrical geometry (Sec. IVIII) . 

In Sec. HV] the humidity dependence of ionic conductance associated with electric field 
confinement within the water layer is calculated and compared with existing experimental 
data. We obtained an excellent agreement between theory and experiment in spite of the 
approximate character of our electrostatic model. A method for experimental verification of 
our model is to measure the non-linear voltage dependent conductance as proposed in Sec. 
IVl We predict that a relatively small electric field ~ lO^V/m, such that FJc ~ ksT is 
needed to sufficiently reduce the ion binding energy to produce a non-linear current-voltage 
dependence. 



II. MODEL 



Consider a partially hydrated DNA molecule shown in Fig. [T] One can distinguish 
dielectric cylindrical layers of DNA, water and air characterized by their dielectric constants 



ein (DNA), ew (water) and eo (air), respectively, and by radii rj„ (DNA) and Tq. The 
radius of DNA and bulk dielectric constants of air and water are known (rj„ ^ 1.3 A in 
an A-conformation, e^, ~ 1 and ew ~ 80). The dielectric constant of DNA is not well 
known because of the obvious difficulty of its direct measurement. Different values for ej„ 
are proposed in different works. A low value for DNA dielectric constant e ~ 2 — 4 was 
suggested in Ref. [28] which is typical of other organic polymers. Further experimental 
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3l|, 



and theoretical studies suggest the dielectric constant of DNA is larger: e ~ 7 
e ~ 12.4 — 20 [3^, e ~ 100 33]. We used different values for and found that the 
intermediate value 4: < ein < 12.4 leads to an agreement of our theory with experimental 
data for humidity dependent DNA conductance. 

The humidity dependent radius of the water layer ro is found using the experimental 
results for the number of adsorbed water molecules per nucleotide, iV^** (see Table 1 in 20|, 
Fig. 1 in Ref. 2l|, 22|]). It is important to note the adsorption measurements ignore the 



zeroth layer o: 



humidity 21 



" water which contains approximately 6 molecules per base pair even at zero 



22 



231 ]. These water molecules are bound to DNA due to electrostatic forces 



and/or hydrogen bonds. The total number of water molecules per base pair can then be 



defined as 



be described 



2N^'^ + 6. Alternatively, the number of adsorbed water molecules can 
3y the Branauer-Emmett-Teller (BET) equation, Ny^^ = 44r/i(l — rh)~^{l + 



351 ]. where rh is relative humidity. In our study, we use the experimental 



relationship because results from the BET equation can deviate from real behavior either if 
the water layer thickness approaches atomic size {rh < 0.1) or is comparable to the distance 
separating different DNA molecules rh —>■ 1. The experimental data has no such limitations. 
The humidity dependence of N\y is summarized in Table [H 

Using the number of water molecules per nucleotide Nw at different relative humidity 
we can calculate the area of the water tube as A]y = Nwivw I dm) and its radius as = 
■\/ Aw/t^ — rf^, where vw = 29. gA^ is the volume per single water molecule and din = ^AA 



is the period of an A-DNA chain. The period of the DNA chain is the average distance 
between base pairs, which varies with DNA conformation. We consider A-DNA because 



according to Ref. 



221 ]. DNA remains in A-conformation at humidity below 93%. The radius 



of the water layer is larger and the binding energy is lower in A-DNA compared to B-DNA, 
so in the case of coexisting phases the formation of conducting ions will be more efficient in 
A-DNA. Correspondingly, we use A-DNA radius rj„ ^ 1.3nm for the inner cylinder in Fig. 
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rh 


0.05 


0.11 


0.23 


0.33 


0.39 


0.49 


0.54 


0.59 


0.65 


0.75 


0.80 


0.84 


0.87 


0.92 


0.93 


Nw, exp. 


8.64 


9.36 


10.6 


12.1 


12.7 


14.5 


15.4 


16.7 


18.2 


22.9 


26.3 


27.1 


32.2 


46.1 


47.5 


Aw, exp. 


1.08 


1.17 


1.33 


1.51 


1.59 


1.8 


1.92 


2.07 


2.26 


2.84 


3.28 


3.37 


4.02 


5.74 


5.91 


To, exp. 


1.43 


1.44 


1.45 


1.47 


1.48 


1.5 


1.52 


1.53 


1.55 


1.61 


1.65 


1.66 


1.72 


1.88 


1.89 


Nw, theor. [35] 


8.38 


9.53 


10.9 


12.0 


12.7 


14.2 


15.2 


16.4 


18.3 


23.3 


28.1 


33.2 


38.6 


60.8 


70.5 


A\Y, theor. 


1.04 


1.19 


1.36 


1.49 


1.58 


1.76 


1.89 


2.04 


2.28 


2.9 


3.5 


4.14 


4.81 


7.57 


8.78 


To, theor. 


1.42 


1.44 


1.46 


1.47 


1.48 


1.5 


1.51 


1.53 


1.55 


1.62 


1.67 


1.73 


1.79 


2.02 


2.12 



TABLE I: Humidity (rh) and water layer parameters N]y, Aw (nm^) and Vq (nm). 



m 

It is not obvious whether the continuous medium approach is apphcable to the ultrathin 
water layers containing few molecular layers (see Table [T]). The electrostatic model is more 
justified if we have at least two or more molecular layers of water, which takes place at relative 
humidities exceeding 0.5. Since our model provides a good description of experimental data 
at all humidities we believe that it can be extended to lower humidities as well. One should 



notice that the recent modeling of dielectric relaxation of cytochrome c oxidase 36|] shows 
that even the ultrasmall hydrophobic cavity around the catalytic center in cytochrome c 
can yet be characterized by the bulk water dielectric constant 80. This gives an additional 
justification for our approach. 



III. ION BINDING ENERGY AND FIELD CONFINEMENT 



Following the previous discussion 



we investigate the ionic conduc- 



tivity associated with H^O^ and 0H~ ions formed in the water self-dissociation, which is 
referred to below as "pure water" conductance. The conductance depends on the density of 
ions, their mobilities, area of the water layer, and the length/ number of DNA molecules. 
We expect that the ion binding energy possesses the most significant humidity dependence 
because ion density depends on this energy exponentially. Only exponential dependence 
can be responsible for producing the six orders of magnitude change in humidity depen- 
dent DNA conductance. The conductance can be approximately expressed through the ion 
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binding energy Ub as jl^ 

Cocexp(-f/B/(2A;BT)). (1) 
The binding potential energy of two ions can be separated into two parts 

Ub = Uw + f/c. (2) 

The first contribution is the direct Coulomb attraction which leads to the binding energy in 
bulk water, Uw- For instance the binding energy of OH^ and ions forming H2O molecule 
is given by Uw 0.48eV [38]. The second contribution to Ub is is due to ion interaction with 
the induced polarization of dielectric environment at interfaces. In the electrostatic model 
the humidity effect is fully related to the interaction with induced surface charges (see Sec. 
IVII|) which create their own contribution to the binding energy Uc absent in the bulk medium 
(ro 00). Changes of humidity modify the thickness of water layer, Tq, and thus affects 
interaction of ions with dielectric polarization. The conductance dependence of humidity 
can be expressed as 

C = Coexp{~Uj{2kBT)), (3) 

where Cq is the preexponential factor describing the conductance in the limit when the 
radius of water layer approaches infinity {Uc 0). The binding energy Uc depends only 
on the ion charges only because it is determined by the induced Coulomb interaction at 
distances larger than its ionic size (see Sec. IVIII) . Therefore humidity dependence should 
be similar for different mechanisms of ion formation including water self-dissociation and 
cation escape from the charged phosphate group (see discussion in Sec. [IV|l . 



If, in accordance with 
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14l |. the conductance is due to 0H~ and H"*" ions formed by 
self-dissociation of H2O, then the factor Cq is determined by the conductivity of a pure bulk 
water {cw ~ 5.5- lO^^Cm/m). The preexponential factor in the "pure water" conductance of 
the junction Eq. ([3]) made of DNA molecules of the length / connecting two electrodes can 
be estimated as Cq = cwAy/N/l so the conductance of the whole system can be expressed 
as 

C = cw^^exp{-Uc/{2kBT)), (4) 

Unfortunately, we did not find any information in experimental literature about the number 
of DNA molecules forming the junction and therefore we are unable to compare theoretical 
prediction for the absolute value of the "pure water" conductance with the experiment. 
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The humidity dependent binding energy Uc can be expressed as the work to bring a 
negative charge, initially located in the direct vicinity of a positive charge, to an infinite 
distance along the cylinder axis. Due to the large difference between the dielectric constant 
of water and those of neighboring layers, the electric field is substantially confined inside the 
water layer L^deed. siuce the uorma. co,„po,>e,>t of e.ectr.c field. E„, on the border of 
two dielectrics satisfies the condition ewEnw = ^mEnm {m = Air or DNA) [37], we can set 
Enw ~ in the limit em/^w 0. In this case, the force lines of electric field (blue arrows 
in Fig. [1]) do not escape from the water tube until the distance between the charge and the 
reference point is greater than the field confinement length Ic < To, thus representing the 
size of the surface charge domain. The electric field, F{r), within the domain restricted by 
the confinement length (r < Ic) can be estimated using the Gauss theorem; the outward fiux 
of an electric field through a closed surface (FcAw) is equal to inkq/e where g = e is the 
charge of the region inside the surface. This yields the following equation 

F(r) = ~ 27ike/{ewAw). (5) 

This distance independent electric field is caused entirely by the surface charges which 
prevent the escape of field lines from the water layer. Therefore, it describes the contribution 
to the binding energy related to DNA hydration that we are looking for. It is important to 
note this field vanishes in the bulk limit Aw — > corresponding to 100% relative humidity, 
where the binding energy is determined solely by direct Coulomb forces. 

The confinement length can be estimated using the continuity of the electric field tan- 
gential component at the boundary between two dielectrics and the field representations 
at confinement length as either the internal field, Fc ~ 27Tke/{ewAw), extrapolated from 
shorter lengths (r < Ic) or the Coulomb field of a point charge, Fc=ke/l1, extrapolated from 
longer lengths (r > Ic). This yields Ic ~ a/ ewAw/i^ir) . Finally, the energy to separate 
charges can be estimated as Uc = eFclc ~ ke^ \/2t: / \/ewAw . Aw ~ Inm^ (see Table H]), we 
can estimate the confinement length as Ic ~ 4nm (accurate calculations result in a value 
one order of magnitude larger, see Fig. [3]). This size is much larger than all characteristic 
molecular sizes, which gives additional justification for the validity of our electrostatic ap- 
proach. Using our estimate for the ion binding energy associated with field confinement, we 
can predict the linear conductance logarithm dependence on the inverse square root of the 
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water layer face area (cf. Eq. (jlj)) 

log(C) = log(Co) - V 



^2ewAwKBT\og{lQi) 
= log{C,)-r^SA/A-^'\ (6) 

where ~ 1 is the proportionahty constant between the estimate and the exact solution 
and Sa = 7.76nm~^. 

Our estimate is very close to exact solutions of Poisson equations in cylindrical geometry 
given in Sec. IVIIt where the direct Coulomb interaction and confinement binding energy are 
separated quite naturally. These solutions are obtained using textbook methods [s^]. The 
solution for binding energy can be expanded over modified Bessel functions of the first and 
second kind, Im and as 

2i-p2 °o f-co Al,K 

Ub{R) = V(2 - 6^0) / dk{¥^lUkR) + ^lK^{kR))- vl/^^ = --; 

Tre^y Jo D 

= {ew - ein){eo - ew)Kra{kR)Im{krin)I'^{krin)Km{kro)K'^{kro) 
-{ew - eo)Im{kR)K^{kro)K'^{kro){einIm{krin)K'^{krin) - twI'm{krin)K^{krin)); 

A^ = (ew - ei„)(eo - ew)Im{kR)Im{krin)I'^{krin)K^{kro)K'^{kro) 
-{ew - ein)Km{kR)Im{krin)I'rn{krin){toIm{kro)K'^{kro) - twI'm{kro)Km{kro))] 
D = {ew - ein){eo - ew)Im{krin)I'm{krin)Km{kro)K'^{kro) 
-{einIm{krin)K'^{krin) - ewI'm{krin)Km{krin)){eoIm{kro)Kl^{kro) - ewl'^{kro)Km{kro)){l) 

Binding energies of ions Eq. ([7]) are calculated using the Matlab software package. This 
energy depends on the position R of the positive charge (dissociating water molecule) with 
respect to the cylinder axis, < R < Tq. The binding energy has a minimum at the 
intermediate value of radius R = R^ (Fig- EI) and tends to infinity at the layer boundaries, a 
characteristic of the point charge model. Since water splitting is most efficient for molecules 
at the binding energy minimum separated from both boundaries by distances exceeding the 
water molecule size, the point charge approach is justified. Thus we used Eq. ([1]) with 
binding energy Ub{R*) to compare our theoretical predictions with the experimental data. 
In the domain of interest R ^ R^, the series of Bessel functions (Eq. (^^) rapidly converges. 
Use of 10 — 20 terms in expansion is sufficient to provide 1% accuracy in calculations. 

The humidity dependence of the binding energy is caused by the related dependence of 
the outer radius Tq given by Table [B This radius is determined by humidity; it tends to zero 
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if relative humidity approaches and goes to infinity when the relative humidity is close to 
1. Using this dependence we calculate the humidity dependence of conductance. 

Note the confinement energy does not depend on the chemical nature of ions formed 
in dissociation, but rather on their charges. Therefore, the formation of conducting ions 
by dissociation of any molecule can be characterized by the same humidity dependence of 
conductance. 




FIG. 2: Confinement energy at different positions of dissociating (water) molecule at relative 
humidity 0.75. 

To describe the field confinement effect, we also consider the distance dependent potential 
energy of two charges. Both charges are located at distance -R* from the polar axis and 
azimuth angle = 0, and separated by distance z. The potential energy interaction related 
to confinement can be expressed similarly to Eq. ([7]) as 

\Jb{R) = V(2-5„o) / dk{^'^I^{kR:)^^'iK^{kR:))cos{kz). (8) 

This expression was calculated using Matlab software for relative humidity 0.5 and the 
results are shown in Fig. [31 It is clear from the graph the confinement related part of 
the interaction dominates at long distances of 25 — 50nm; therefore, the confinement length 
exceeds lOnm. Since the associated potential barrier length for charge separation is very 
long, this barrier is extremely sensitive to the external electric field so strong non linearity 
in current- voltage dependence can be expected. This effect is discussed later in |Vl 
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FIG. 3: Point charge potential dependence on the distance; rh=0.5. Induced field clearly exceeds 
the Coulomb field at distances exceeding few A. Electric field Fc = 2.5 • 10^ V/m reduces the 
potential barrier for water self-dissociation by the thermal energy /c^T ~ 0.026eV while the field 
10''V/m destroys confinement effect. 

The dependence of conductance on the water layer area (Eq. ([T])) for different DNA 
dielectric constants is shown in Fig. |H The calculated behavior agrees with the qualitative 
prediction (Eq. ([6])). Moreover, the slopes of the linear conductance dependence on A^' 
also agree with Eq. (5) Eq. (Table HTj) . which proves the significance of an electric field 
confinement. 




FIG. 4: Dependence of conductance on the area of water layer. 
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IV. COMPARISON OF THEORY WITH EXPERIMENTAL DATA 



Combining calculations of conductance dependence on water layer area (Fig. H]) and 
experimental data for area dependence on relative humidity (Table [T]), one can describe the 
humidity dependence of conductance. This dependence is shown in Fig. O for different 
dielectric constants of DNA (ei„) together with the available experimental data. We were 
not able to predict the pre-exponential factor; therefore, all conductance values are given in 
arbitrary units. 

It follows from Fig. [5] that the calculated conductance shows exponential dependence on 
humidity as in experiments; moreover its behavior agrees quantitatively with majority of the 
experimental data for 4 < ej„ < 12. The theoretical and experimental slopes of exponential 
conductance dependence on humidity are similar as shown in Tab le [ITl We can see if the 



DNA dielectric constant is set to 7 as suggested in Refs. 29, |31| . then we obtain an 
exponential dependence slope dlnC/drh ^ 0.58. This value is identical to the one observed 
in several studies including poly-A poly-T DNA [Qi], single and double stranded herring 
DNA Q, and A-DNA_|lJ In other experiments with poly-G poly-C DNA 9] and DNA 
linked gold particles [l2|, larger slopes (0.7 and 0.9) have been reported. Much stronger 
dependence has been observed for DNA modified with acridine orange |lo| which is most 
likely due to its strong chemical difference from standard DNA. Our electrostatic model can 
fit all of this data (Table [IT|) . 

Note that the slopes can be sensitive to fluctuations in the bundle structure supporting 
water dissociation. The obvious example of such fluctuation is the overlap of two adjacent 
DNA molecules or the water layers formed around them. In this structure, the water layer 
face area Aw will be approximately doubled (cf. Fig. [1]). This doubling reduces the 
humidity dependent part of binding energy by a factor of -\/2 (([6])), which also reduces the 
slope of the conductance exponential dependence by the same factor. The overall effect 
of fluctuations depends on their probabilities which is not known to us, but can vary in 
different experiments. The fluctuations will lead to the weakening of humidity dependence 
and the experimental data characterized by the weakest dependence could have the strongest 
fluctuation effect. 

In addition to the humidity dependence, other interesting properties of DNA conduc- 
tance at various humidity were investigated, including the sequence and temperature de- 
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pendencies [9|, |12| . and comparison of humidity dependent conductance of single and double 
stranded DNA [13*]. The experimental temperature dependence of DNA conductance has 
been well approximated by the Arrhenius law. In Ref. [ol, the temperature dependence 
is reported at relative humidity 0.5 and the conductivity activation energy has been esti- 
mated as 0.49eV for a poly-G poly-C sample and 0.6eV for poly-A poly-T sample. We 
can estimate the ion transport activation energy, Ea, for "pure water conductance" using 
Eqs. ([1]), ([2]) and our calculations. For DNA dielectric constants = 2,4,7,12.4,20 we 
get, Ea = 0.62, 0.6, 0.58, 0.55, 0.52eV, respectively. All estimates agree reasonably well with 
experimental values. 



In Ref. 



12l |. the activation energy was studied for large relative humidities, 0.7 and 0.85. 



It was found to be less than that of pure water conductance which is 0.24eV. Therefore, 
the formation of ions in this system is most likely associated with a process other than the 
self-dissociation of water molecules. Yet this process seems to also be sensitive to the field 
confinement and humidity in good agreement with our theory. Perhaps the process is the 
breaking of a bond between the ion and phosphate group which has smaller bulk binding 
energy than that of the water molecule. 

The humidity dependent conductances in single and double stranded herring DNA [Ic 
are found similar to each other. This is not what is expected in our model. According to 



Ref. 



2l| , the number of adsorbed water molecules per nucleotide is approximately the same 



for both single and double stranded DNA. Accordingly, the total number of water molecules 
per unit length and area of the water layer. Aw, is smaller for single stranded DNA by a 

1/2 

factor of 2. Since the confinement energy Eq. ([6]) is inversely proportional to A^r the 
binding energy associated with the field confinement should increase by the same factor of 
and so does the slope of the humidity dependence. This expectation confiicts with the 
experimental observation of nearly identical slopes for single and double stranded DNA. 
How can this discrepancy be explained? It is possible the conductance was measured in an 
essentially non-linear regime in Ref. 13] (current voltage dependence has been reported at 
bias voltage 4V where nonlinearity is obviously significant). Our theory is only applicable 
to the linear regime. One should notice the four-contact measurements in Refs. 
are performed in the linear regime, while nonlinearity is clearly seen in the current-voltage 



characteristics observed in Refs. 



12| . Perhaps the small value of slope and the identical 



behavior of single and double stranded DNA conductance is the consequence of non-linear 
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2 


4 


7 


12.4 


20 


dloq(C) ,1 

V = „ ," -1/2. theory 


1.49 


1.39 


1.25 


1.06 


0.87 


Slope ^^^7^, theory 


0.7 


0.64 


0.58 


0.49 


0.4 


Experiment 


[9] 




[12j 


[13j 


M 


Details 


GC 


AT 


DNA linked gold 


herring DNA 


A-DNA 


Slope (rh), experiment 


0.76 


0.58 


0.9 


0.58 


0.58 



TABLE II: Conductance dependence on the area of water layer and humidity 

regime of measurements. Also, the two contact method was used in both Refs. Q,Q. The 
conductance measurements made using this method are sensitive to the junction resistance, 
which is excluded in four contact measurements jol. 



O 
o 



-^2■ 




0.0 



0.5 

Relative humidity 



1.0 



FIG. 5: Experimental data and theoretical predictions for humidity dependent conductance. 



V. DISCUSSION AND CONCLUSIONS 

Thus we suggested the solution of the long-standing problem of the humidity dependence 
of DNA conductance. The dramatic increase in DNA conductance of six orders of magnitude 
with increasing humidity is interpreted using the change of ion binding energy induced by 
the electric field confinement within the water layer. This confinement is caused by the 
large value of the water dielectric constant compared to the environment. Consequently the 
binding energy of conducting ions is very sensitive to the thickness of water layer determined 
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by humidity. Increasing the humidity results in a growing water layer, which reduces the field 
confinement. Consequently, attraction force between two oppositely charged ions decreases 
and activation energy of ionic transport decreases accordingly leading to a dramatic increase 
of ionic density and conductivity. This behavior is opposite to the one expected for electronic 
transport where activation energy is expected to grow with the amount of water surrounding 
DNA 0. 

Yet depending on particular experiments, there are problems in interpretation of con- 
ductance mechanisms due to strong variations in conductance temperature dependence [3] 
and the absence of denaturation effect of conductance [l^- Therefore, it would be inter- 
esting to perform direct experimental verification of field confinement. We believe that the 
measurements of a non linear current-voltage characteristic can be used as such verifica- 
tion. Strong nonlinearity of current-voltage characteristics of DNA is expected because the 
confinement length is much longer than characteristic molecular sizes (Fig. [3]). Indeed the 
external DC electric field F should reduce ion binding energy and support ion dissociation 
in the relatively large domain of field confinement restricted by confinement length Ic (see 
Fig. [3]). Assuming this DC field is much smaller than the confinement field Eq. ([5]) one 
can describe its effect as a perturbation. Then we can estimate the reduction of ion binding 
energy as 5Uc = —Flc. The nonlinear conductance can be approximately represented as 

C(F)=C(0)e.p(^). (9) 

The nonlinearity should become strong at Fc ~ 2kBT/lc. According to our solution for 
ion binding energy distance dependence (Eq. ([8]), we can estimate Fq ~ 8 x lO^V/m for 
relative humidity 0.5. The threshold field is similar for other values of humidity. The 
measurements of this dependence (Eq. (Q) are very interesting because they result in 
the direct determination of the confinement length. One should notice the four contact 
measurements in Ref. jol were made at DC fields less than 8.3 • lO^V/m. The nonlinear 
threshold for those measurements can be higher than our prediction because DNA molecules 
were aligned under some angle to the field. In Ref. |l2|, the maximum electric field was 
near 3 ■ lO^V/m. We suggest an increase of the field by one order of magnitude which should 
lead to remarkable nonlinearity according to our model. The application of a field exceeding 
20Fc should remove the confinement effect as illustrated in Fig. [31 The exponential field 
dependence should become weaker when the DC field approaches the confinement field Fc. 
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Thus a non-linear current voltage characteristic can be used to extract both confinement 
field and confinement length that can be compared to our electrostatic model. 

n 

Note added in proof. While our paper was under review, another paper [39] considering 
humidity dependent conductivity of DNA appeared in press. This paper suggests an expla- 
nation of experimental data using two-dimensional field confinement in planar water layers 
surrounding DNA films. Their model is both quantitatively and qualitatively different than 
our model which uses one-dimensional field confinement. We believe that further experimen- 
tal investigation of humidity dependent conductivity, particularly in the nonlinear regime, 
will help choose which explanation is more relevant. 
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VII. APPENDIX 



Here we derive the exact solution for the potential of a point charge located in the water 
layer between DNA and air Fig. [1] and express the ion binding energy using that potential. 
The charge e is placed in water and its cylindrical coordinates are chosen as ro = (p, 0,0), 
fin<P< fo- The charge potential satisfies the Laplace equation 

A$ = 47rA;-(5(r-ro), (10) 

and the standard boundary conditions of continuity of potential in the border separating 
different media and discontinuity of its normal derivatives 

,^^ = ^2^ (11) 

drii dn2 ' 

Using the cylindrical symmetry of the problem the solution can be expanded into the series 
of modified Bessel functions following the textbook [s?!. Inside the DNA (r = (p, z, ip)) one 
should use modified Bessel functions of the first kind having no singularities at p = 0. The 
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solution at the point < r < rj„ can be expressed as 

2ke 



poo 

V(2-5„o) / dk¥;^{k)lUkp)cos{kz)cos{m^). (12) 

, n ^0 



where Smo is the Kronecker symbol and functions unknowns to be determined. 

The solution in the air < r can be expressed similarly with the only difference that 
one should use modified Bessel functions of the second kind approaching zero when their 
argument tends to infinity 

2kp 

^outir) = V(2 - 5mo) / dk^':;:\k)K^{kp)cosikz)cosimip). (13) 

The solution inside the water layer can be expressed as the superposition of the field of 
point charge satisfying Eq. (ITOl) and the general solution of Laplace equation using modified 
Bessel functions of both first and second kind (cf. Ref. [s^]) 

?kp ^ 
<|.^(r) = ^5^(2-M 

POO 

X dk {Km{kpo)Im{kp)6{po - p) + Kjn{kp)lm{kpo)0{p - Po)) cos{kz)cos{mip) 
Jo 

2ke 



+ 



y"(2-5,r^o) / dk{¥^{k)Im{kp) + '^^{k)K„,{kp))cos{kz)cos{m^). (14) 

™_n Jo 



The first term represents the Coulomb field of the point charge in the dielectric medium, 
which is independent of the system geometry and describes the case of an infinite homo- 
geneous sample. The second term is the surface charge contribution which describes the 
geometry effect and correspondingly, the humidity dependence of interest. This contribu- 
tion to the binding energy is determined by the potential difference between the bound 
state, ^^^(ro) (the position of the second charge coincides with the position of the first), 
and the state where the second charge is infinitely far from the first, $vy(oo) = 0. Thus the 
confinement induced binding energy for ions having charges e and — e can be expressed as 

Uc = V(2 - 6mo) / dk {^l{k)lUkp) + ^l{k)K^{kp)) . (15) 

^0 Jo 

Unknown expansion coefficients ^'^^ can be determined using boundary conditions of poten- 
tial continuity and electric field discontinuity Eq. (iTTl) at the boundaries. These expressions 



16 



along with the binding energy Eq. ([7]) are given within the main body of the present text. 
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